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Abstract 

We study condensation in several particle systems related to the 
inclusion process. For an asymmetric one-dimensional version with 
closed boundary conditions and drift to the right, we show that all 
but a finite number of particles condense on the right-most site. This 
is extended to a general result for independent random variables with 
different tails, where condensation occurs for the index (site) with 
the heaviest tail, generalizing also previous results for zero-range pro- 
cesses. For inclusion processes with homogeneous stationary measures 
we establish condensation in the limit of vanishing diffusion strength 
in the dynamics, and give several details about how the limit is ap- 
proached for finite and infinite systems. Finally, we consider a con- 
tinuous model dual to the inclusion process, the so-called Brownian 
energy process, and prove similar condensation results. 

Keywords: inclusion process, condensation, Brownian energy pro- 
cess, zero-range process. 

1 Introduction 

In PQ, 0, an interacting particle system was introduced, where particles 
perform random walks and interact by "inclusion", i.e., every particle at site 
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i can attract particles from a site j to its site at rate p(i,j) = p(j,i). This 
particle system, the so-called symmetric inclusion process (SIP), is "exactly 
solvable" by self-duality, and its ergodic stationary measures are products of 
discrete gamma distributions, indexed by the density. The inclusion process 
also turns out to be dual to a system of interacting diffusions, the so-called 
Brownian energy process (BEP). More details on duality, self-duality, and the 
precise relations between SIP and BEP can be found in [Tj. In the present 
paper we only need the explicit form of the stationary measures of these 
models. 

We prove existence of stationary product measures for inclusion processes 
under rather general conditions, in analogy to classical results for exclusion 
processes [TB] . We introduce asymmetric versions of the SIP and the BEP, 
for simplicity focusing on a one-dimensional context with N sites and closed 
boundary conditions. In this case both models have spatially inhomogeneous 
product measures as reversible measures (to be compared with the blocking 
measure of the asymmetric exclusion process). Conditioning on K particles 
in the system (resp. total energy E), we prove that that in the limit K — > oo 
"almost all" the particles (resp. all the energy) are concentrated on a single 
site, where the marginal of the reversible measure has the heaviest tail. The 
other sites contain a finite number of particles (resp. finite amount of energy). 

We further study condensation in inclusion processes with spatially ho- 
mogeneous stationary measures, with the SIP as the main example. The 
strength of the diffusive part of the dynamics in comparison to the attrac- 
tion is controlled by a system parameter m > 0. For fixed particle density 
p we study the limit m — > where attraction dominates, and show that the 
single-site marginals converge to Dirac measures concentrated on zero mass. 
This corresponds to the fact that a typical configuration consists of rare piles 
of typical size 2p/m separated by empty sites. The distribution of pile sizes 
approaches a power law with exponent —1 and becomes degenerate in the 
limit m — > 0. This leads to a breakdown of the usual law of large numbers 
which we illustrate in detail. 

Our results for the asymmetric case also cover condensation phenomena 
in zero-range processes, which have attracted a lot of recent research interest 
[SI E|. For inhomogeneous systems, these have been studied before mainly 
in the context of a quenched disorder in the jump rates, which have to be 
non- decreasing functions of the number of particles [HI El HD]- For such 
systems, the use of coupling techniques allowed in special cases to also obtain 
results on the dynamics of condensation. In contrast, our results cover only 
the stationary behaviour but apply to a much larger class of jump rates with 
essentially no restriction. The widely studied condensation in spatially ho- 
mogeneous zero-range processes [HI [T2J H31 El [15] has a somewhat different 
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origin than our homogeneous results for the SIP. This is discussed in detail 
at the end of Section 4.2. 

In the next section we describe the inclusion process and its stationary 
measures. In Sections 3 and 4 we study condensation in the asymmetric and 
spatially homogeneous case, and discuss extensions and relations to zero- 
range processes. In Section 5 we introduce the asymmetric Brownian energy 
process and discuss condensation in an example of a system with continuous 
state space. 

2 Inclusion processes 

The inclusion process on a general discrete set A has state space Q = N A and 
we denote a configuration by rj — (rji : i e A) where rji is interpreted as the 
number of particles at site i G A. The dynamics is defined by the generator 
defined on the core of local functions / : f2 — > R: 



where r\ % ^ is the configuration obtained from 77 by removing a particle from 
site i and putting it to j. The p(i,j) > are jump rates of an irreducible 
random walk on A with p(i, i) = 0, and the parameter m > determines the 
rate of diffusion of the particles as compared to the aggregation part given 
by the product rjiTjj. We also assume the p(i,j) to be uniformly bounded and 
of finite range, i.e. there exist C, R > such that 

sup p(i,j) < C and | { j e A : p(i,j) > 0} | < R for alH e A . (2) 



This ensures that the dynamics is well defined even on infinite lattices (for a 
large class of 'reasonable' initial conditions) and contains all generic examples 
we are interested in, such as nearest-neighbour hopping on regular lattices. 
If the p(i,j) are symmetric the inclusion process is also called symmetric 
(SIP), otherwise asymmetric (ASIP). 

2.1 Stationary product measures 

For > and Aj > 0, i G A, define the product probability measure 



Lm = J2 p&fi* (y + »&) - /fa)) > 



(1) 



i,jeA 



ijeA 



(3) 
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where the marginals v % are probability measures on N given by 

#H^(f^T^ (4) 
with the normalizing constant 

^m=f:^ T -^^=(i-^)- m/2 . (5) 

n=0 ' LL \ 2 ) 

The parameter <fi > is called fugacity and controls the particle density, 
which is invariant under the time evolution. 

Theorem 2.1. For all <fi < <p c := (sup igA Aj) 1 , is a stationary measure 
for the inclusion process with generator (T7]) ; provided that one of the following 
conditions holds: 

a) Thep(i,j) are doubly stochastic modulo a constant, i.e. 

^2(p(hj)-p(j,k))=0 foralli,keA, (6) 
jeA 

and Xi = 1 for all i G A. 

b) The Xi are reversible w.r.t. thep(i,j), i.e. 

\p{hj) = AiP(j,z) for all i,j G A , (7) 

and in that case is also a reversible measure. 

This is in direct analogy with well-known results for stationary measures 
for exclusion processes (see e.g. [16], Thm VIII. 2.1). In both cases, the Aj 
are special harmonic functions solving 

( Xi P& ~ X iPti> *)) = f or all z G A , (8) 

jeA 

i.e. they provide a (not necessarily normalized) stationary distribution for 
the underlying random walk of a single particle. For the above product mea- 
sures to be stationary, the p{i,j) have to be such that they admit a constant 
solution (first case) or a detailed balance solution (second case). It is not 
clear at this point whether these conditions are really necessary for the ex- 
istence of stationary product measures in general. Note also that on infinite 
lattices 4> c = is possible. But for finite A (which we mainly focus on in 
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this paper), Theorem 12.11 guarantees the existence of a family of stationary 
measures. 



Proof. We have to show for expected values w.r.t. that 

MW) = E E P(»*. J)* (y + - = (9) 

?jGO ij'eA 

for all local functions /. For fixed i,j we get after a change of variable 

rten 

The form (J4]) of the marginals implies that for alH £ A and > 

+ _ m + 2fc 

^2(fc + l) ^ • 

Thus we get for each fixed pair z, j £ A 



^ n+1 



) 4(k-l) (n + 1) (y + A; - l) = i^(fe) fc + n 



A, 



for all n > and A; > 1. It is easy to check that boundary terms in the sums 
vanish consistently, and we do not consider them in the following. Plugging 
this into flU]) we get 



(10) 



and exchanging the summation variables i <H- j in the first part of the sum 
leads to 

MLf) = Y,f(vH(v) E Vi{m/ l + Vi) (pU,i)*i-p(i,m ■ (11) 

This clearly vanishes under the reversibility condition ([7]) which implies sta- 
tionarity under assumption b). In analogy to ffTOl we can derive 



v<t>{gLf) = ^2f(v)Mv) 



E (y+^) yaW) ~ Vi ^ 



and after using (j7|) and the exchange of summation variables this implies 

= Mf L a) , 

so is also reversible. 

Assuming a), the A, and Xj in (jlip cancel, and we write the linear (diffusive) 
part as 

2-5^(p0'> 1 ') -Pihj)) = > 
which vanishes due to <$5§. For the quadratic aggregation part we get 

Z (p0'> *) _ Pfcj)) = ViVj (pU, *) - P(j\ *)) = ' 
i,ieA tjeA 

by another exchange of the summation variables in the second part, using 
that r]ii]j is symmetric under % •<->■ j. □ 



2.2 Canonical measures for finite systems 



Consider a finite lattice A at of size N with corresponding state space Qn = 
N An . Starting with a fixed number of K particles, the inclusion process with 
generator as given in (TjQ) is an irreducible continuous-time Markov chain 
on the finite set 



N 



A K = < 7] e Q N : Y Vi 



K 



(12) 



i=i 



and has a unique stationary measure, which we denote by fi . 

By conservation of the number of particles, the conditional measure 



dr] 



N 



K 



8=1 



is also invariant. Indeed for / : Qn — > K we have 

f L N f(r))l AK M d v) 



K 



Jf( V ) (L* N (1 Ak )) (^(dv) 



v<t>(A 



K 



(13) 
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since it is easy to see that with the generator L N also its adjoint L* N con- 
serves the number of particles. In the case of reversible measures u^, is 
self-adjoint and there is nothing to check. By uniqueness of the stationary 
measure, we thus have 

^(. \A K ) = /2 K (14) 

for all (ft < 4> c and K £ N. So the conditioned product measures are actu- 
ally independent of 0, and this connection provides an explicit form for the 
canonical measures \i K . 

3 Condensation in the ASIP 

A generic situation where Theorem 12. II gives rise to spatially inhomogeneous 
reversible measures is a one- dimensional lattice = {1, . . . , N} with an 
underlying asymmetric nearest-neighbour walk. We consider the ASIP with 
generator 

N-l 

LnM = E^(Y + r?m ) (/(r/M+1) " /(r?)) 

N-l 
i=l 

where p > q > 0. In this case Aj = {p/q) 1 fulfills condition (JTJ) in Theorem 
12.11 We will now proceed towards showing that in the limit K — > oo, under 
the canonical measure p, K , the typical situation will be that all but a finite 
number of particles condenses at the right site i = N, whereas the other sites 
contain a number of particles distributed according to v^ c . 

At first sight one could be tempted to think that this is just a consequence 
of the asymmetry: particles are pushed to the right. This is, however, not 
the case. If we consider independent random walkers, moving at rate p to 
the right and q to the left, then the reversible profile measures are Poissonian 
and given by <8>iLiV\(dr)i) with 

with a normalizing constant Zi(<p) = e^/vY which is now finite for all values 
of (j). As a consequence, no condensation happens: if we condition on having 
K particles, and let K tend to infinity, all sites will carry a diverging number 
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of particles. The condensation phenomenon is thus a combination of the 
asymmetry, together with the attractive interaction between the particles in 
the inclusion process. Indeed, it is the interaction which is responsible for 
the existence of a finite critical </> c . 



3.1 Condensation 

Before we formulate the main result of this Section, we recall the marginals 
v\ for the ASIP 

^n) = ^—r^{n) , (16) 

where we have now 

\ = (p/qY (17) 

and write 

<n) = ^-tfi . (18) 
nil (n) 

In the present case Wi does not dependend on i, but in generalizations ex- 
plained below we will allow explicit dependence on i. The weights Wi(n) have 
the asymptotic behavior 

Wi(n) ~ n 2 (19) 

where means that a n /b n converges to a strictly positive constant. 

We remind that the normalizing constants are 

-m/2 

Zl (<p) = (i - \ l( pr m/2 = (i- ( ?.)'<!>) . (20) 



Therefore, in the context of Theorem 12.11 we have Ai < A2 < ... < \n, 
(fi c = 1/Atv, Zi((f> c ) < 00 for all 1 < i < N — 1, and z N (<fi) < 00 for all <fi < <p c . 
We then have the following result. 

Theorem 3.1. a) In the limit K — > 00, rji, . . . , are asymptotically 
independent and converge in distribution to the critical product mea- 
sure, i.e. for all m, . . . , njv-i £ N 

H K (771 = m,.. . ,r)N-i = rin-x) ^J c (ni) ■ ■ •^ c " 1 (^at-i) (21) 

where (f> c = I/Xn = (q/p) N - 

b) In the limit K — > 00, the right edge contains "almost all" particles, i.e., 
for all 5 G (0, 1) 

ti K {riN< (1-5)^)^0, (22) 
and we have a strong law of large numbers, r/^/K — > 1 a.s. . 
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Proof. We use that fi = v${. \A K ) and write for A'C A 



N 



Z(A',K)= Yl Jlwiim))?. (23) 

K,i€A':£ igA , ni=K} i£A' 



We then have 

H K (rjx = m, . . . , 77^-1 = n„_i) 



N-l \ N-l 



Z(A N ,K) 
We first prove that 



w N 



i=l / i=l 



^-E n A ^ 'n^)^- ( 24 ) 



iim §^=n*. (25) 

To see this, we choose an appropriate order of summation, 

K K—n\ K-(ni+...n N - 2 ) 

z(a n ,k) = y: E ••• E 

ni=0 n2=0 ri]v_i=0 

'JV-1 \ / N-l 

n w jM x ? ) w n [ k - e ^ ) a 

oo oo 



N 



a^wE- E 

ni=0 n]v_i=0 

* K (ni, . . . , njv-i) JJ ^K) (^J (26) 



with 



(a- Ejli 1 n j) 



• •) - Ki<K ■ ■ ■ K N -i<K-ni-..-n N - 2 ■ ( 27 ) 

W N {K ) 

We see from ffl8l) that ^/^ < 1. Therefore, by dominated convergence, using 
that (fi c = A^ 1 and 



oo oo / A ■ \ n 

kyc) = E M n ) K = E I A / < ^ 

n=0 n=0 V ' l X ' 



we obtain (J25J). Combining 021]) and (J2"5|) with the fact that 

lim WN ^ K ~ n) = l f or all n G N (2* 

K^oo W N (K) 

(which follows immediately from ()19p ). yields item a) of Theorem 13.11 
To prove item b), we start with 

En<(i-6)K w N (n)X%Z(A N \ {N}, K-n) 



H K (r) N <{l-8)K) 



Z(A N ,K) 

and estimate, for n < (1 — 5)K and a small enough e' > to be chosen below: 
Z(A N \{N},K-n) < (X N ^(1 + e')f- n 



K-n K-n-{n 1 +...+n N ^ 3 ) / N-2 

■ ■ ■ Yl wn -i [ k - n - H 

ni=0 n JV „ 2 =0 V j=l 

'N-2 

] I w A n j) 



A, 



W N -i 



s X N ^(l + e') 

< C (X N ^(1 + e')f- n (1 + e) K . (29) 
Here we have used that (cf. (TP9"j) ) 

N-2 \ 

K-n-J2 n i) < c ( 1 + e ) K (3°) 

for some e > to be chosen below, and the fact that the remaining sums in 
the RHS of fT29|) converge to a finite value as K — )■ oo. By f )25|) Z(A^,K) is 
bounded below by CX^wn^K) for K large enough. This then gives 

since for the summation indices K—n > 5K. Choosing e, e' > small enough 
such that 

< ^ < q < 1 

and using that WN< $^ < 1 , we obtain 

& w N (K) — ' 

l* K (VN < (1 - S)K) < C"q SK . (31) 

Choosing 5 = 5k = 1/vK — > 0, we get a summable bound on the right-hand 
side. Since by definition r] N < K a.s. under the measure \i K , this implies al- 
most sure convergence and the strong law i]n/K — > 1 by Borel-Cantelli. □ 
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3.2 Generalizations 



Notice that in the proof of Theorem 13.11 we did not use the specific form of 
Wi and Aj. Therefore, the same proof shows a condensation phenomenon for 
a general family of independent random variables rjx, . . . , T]n with 

under the following hypotheses on the Wi, Xf. 

a) The Aj satisfy 

N—l 

Xn > maxAj , (32) 

i=l 

b) the weights Wi(n) are subexponential in the following sense 

lim = 1 (33) 

for all 1< i < N. 



From fl33l) . ( 128)1 follows directly, and it further implies that for all a > there 
exists C a > such that for all n > 

C-V " < wiin) < C a e an . 

From this bound we conclude that for all > 0, there exists (Tg > such 
that for all n, I > 0, 

r -l -/3(n+Z) < < /7 a J?(n+Q 

This is all we need in the dominated convergence argument to bound ^ k of 
(12T)1 . and to conclude (1317)) . ( 13T)) . Therefore, under the assumptions a), b) we 
conclude the statement of Theorem 13.11 with 



jV 



Example: Zero-range processes 

Consider a general zero-range process on CIn = N^ 1 '""'^ with generator 

L N f(v) = J2p(hj)9i(Vi) W J ) ~ M) , (34) 

i,0 
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where p(i,j) are rates of an irreducible continuous-time random walk on 
{1, . . . , iV} and where : N — >■ [0, oo) with gi(ri) = if and only if n — 0. 
Moreover, we assume for the moment that gi(n) — > 7j G (0, oo) as n — ?• oo 



for all i = 1, . . . , JV. 



By irreducibility of p(i,j), up to multiplicative constants there exists a 
unique function «:{!,..., iV} — » (0, oo) such that for all z, 



("iKM) ~ K iP(hj)) = • (35) 



i=i 



Under these conditions it is well known [T?] that the zero-range process has 
stationary product measures with marginals 

1 S n K n 

* B) = ^n&w • (36) 

which are of the form ([TED with 



T' 1 

= Tin - 7m and A< = — . 
So in order to apply the general result, we need 



UN JV-l /Cj 

— > max — 



7jv i=1 7i 

and the subexponentiality condition on tfj follows since 

Wi(n+1) 7i 
lim — — — = hm — = 1 . 

n->oo W%[n) n^oo g i [n + 1 ) 

Remark 3.1. 1. The case 7^ = 00 for some i ^ N can be included as 
well. In that case, Zi(<fi) < 00 in ([36D for all <f> > 0, in particular for 
<p = C = 1/ \jy ■ Therefore the result of Theorem \3.1\ still holds. 

2. If there are more sites i such that \ = X^, then a) of Theorem \ 3.1\ 
holds for all i where Aj < Xn- Item b) becomes that all but a finite 
amount of mass is concentrated on the sites where Aj = Ajv- 

Note that we make no assumptions on the jump rates of the zero-range 
process except a regular limiting behaviour, in particular there are no mono- 
tonicity assumptions. The latter have been in place in previous work on 
inhomogeneous zero-range condensation where the gt are non- decreasing [?[ 
El QUI E] , which made it possible to make much stronger statements including 
also the time evolution of the condensation. In that sense Theorem 13.11 is a 
generalization of previous results regarding only the stationary distribution. 
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4 Condensation in homogeneous 
inclusion processes 

In this section we study condensation in spatially homogeneous systems. 
There are two natural situations where Theorem 12. II leads to spatially homo- 
geneous product measures v^. If the p(i,j) are symmetric, i.e. p(i,j) = p(j, i) 
for all i, j G A then the reversibility condition ([7]) is fulfilled by taking a con- 
stant Aj = 1 for all i £ A independent of the geometry of the lattice. The 
same solution holds for translation invariant, asymmetric processes according 
to condition where 

p(hj) — <?(j — *) f° r some q : A — > [0, oo) with bounded support . 

In the second case the lattice also has to be translation invariant, such as 
A = Z d or a finite subset with periodic boundary conditions. The measures 
v$ are then not reversible and the system can support a non-zero stationary 
current of the form 



</(p) = p(y+p)E^) • 



fceA 



4.1 Stationary measures 

In both cases discussed above the inclusion process has a family of homoge- 
neous stationary product measures with marginals 

v >'^w/Sw' (37) 

and the normalizing constant 

Z {4>) = y r T{ 7t n ) = a - <py m/2 . m 

n=0 \ 2 / 

The measures are well defined for all positive <fi < <p c — 1, and the average 
number of particles per site is given by 

Pm((f>) = 0<9 log z((f>) = — ^—^ ■ (39) 

Inverting this relation (f) m (p) = m /2 +p allows us - with a slight abuse of 
notation - to index the measures by the density, 

„ r . (B)= ^w^vem (40) 



z((f) m (p)) \m/2 + pJ n!r(f) 
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We also replace the superscript since the marginals are site-independent and 
we want to stress the dependence on the parameter m. Since the density can 
take all values between and oo, we see that for fixed m > the attraction 
between the particles is not strong enough and the inclusion process does 
not exhibit condensation. However, if we increase the relative strength of 
the attractive part in the generator ([T]) by taking m smaller and smaller at 
a fixed density p, a condensation phenomenon occurs in the limit m — > 0. 



Theorem 4.1. As m — > 0, we have for all p > 

m/2 

!/("»> (0) = 



— \ (l + o(l))-U and 

(\ m/2 / \ n 

i) «"" 2 -(l + O (D)-0 (4!) 

for n > 1, which implies 



m n 
Proof. By direct computation we get that 



"£ m, (w)->r- ( 42 ) 



(\ —m/2 / r, \ m/2 

= - (l + o(l)) ^1 (43) 
m/2 + p J \m J 

as m — > 0, which directly implies the statement for iSp m \o) = 1/z (<j) m (p)). 
For every fixed n > 1 we have 

= I--) (l + o(l)) -^1, 

and using r(x) ~ ■= as x — >■ and (|T9l) we obtain 
r (f + n) m 



n 



ir (f ) ^" 2 



n m/2-l(l + (l)) ^0 , 



which implies the second statement. The limit in ( l4"2"j) follows immediately 
from the asymptotic behaviour. □ 

Therefore, for small diffusion rate m sites are either empty with very 
high probability, or contain a large number of particles to match the fixed 
expected value p > 0. From theorem 14. II we infer the following leading-order 
behaviour for small fixed m, 



p [n) - 2 \ (1 - ^)"n m / 2 " 1 , n > 2p/m 
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n 



Figure 1: The scaled marginal {2/m)u p m for p = 1 and for several values of 
m (full colored lines). The asymptotic behaviour as given in ( T4"4"|) is indicated 
by dotted lines. 

where we have used 



with the notation a m ~ b m if a m /b m — > 1 as m — > 0. 

So the marginals show an approximate power law decay with exponent 
— 1, until an exponential cut-off sets in at the scale n ~ 2p/m. This is 
illustrated in Figure [TJ where we see that the asymptotic behaviour for large 
n fits very well also for smaller values of n. Despite the small prefactor m/2 
the density p > is realized by the asymptotic heavy-tail behaviour, and for 
each m > the distribution is normalized due to the cut-off. Conditioned 
on a site being non-empty, its distribution is given by Up(n)/ (l — Up m \o)). 
Using that to leading order 



Like in (H2l) . convergence is clearly non-uniform due to the cut-off, and the 
limit is not a probability distribution. The interpretation of this result in 





we get with for the conditional distributions 




(45) 
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terms of condensation depends on the geometry and is different for finite and 
infinite lattices A, as discussed below. 



4.2 Finite systems 

For finite lattices one can condition on the total number of particles in the 
system, defining the canonical measures as in Section 2.2. The basic features 
of this approach can already be understood on a system with two sites and 
A = {1,2}. Let 771, 772 be two random variables each distributed as Vp and 
consider their joint distribution /i^ conditioned on their sum being equal to 
K e N, i.e. 

f i«:=v p (.\r ]l + r ]2 = K) . (46) 

For each K e N and m > the inclusion process is irreducible and is the 
unique stationary measure (cf. (fT3l) ). A first observation is that, as before, 
fj,^ does in fact not depend on p since due to cancellation 



K I \ tin 1+ n 2 ,K (n^V^ (n 2 ) 

"-<* = »"* = »"> = 2-^(0^(^-0 

5 ni+n2)K T{m/2 + wi)r(m/2 + n 2 )/ (ni!n 2 !) 
" E£o r ( m / 2 + l)V{m/2 + K- l)/(l\(K - 1)1) 

Proposition 4.1. In the limit m — > we have for all K > 



(47) 



/^^(<W) + <W)K (48) 
i.e. all particles concentrate on one of the sites with equal probability. 
Proof. With 772 = K — r\\ we have 



K 1 K v r (f + + K-n)/(n\(K-n). J 

uz. (m = n, r? 2 = K - n) = — | . (49) 

E£o r (f +0r(f +K-l)/(l\(K - /)!) 



In the normalizing sum, as m — )• 0, the two terms for / = 0, K diverge like 
F(m/2)/K, whereas the rest of the sum converges. Also the term in the 
numerator of ^(771 = n) diverges like T{m/2)/K if n = or K and is finite 
otherwise. This implies the result. □ 



The interpretation is that as m — > aggregation dominates more and 
more over diffusion and the particles tend to cluster on one of the lattice 
sites. The onset of condensation for small m can be well illustrated in the 
limit of infinitely many particles. 
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0.0 0.2 0.4 0.6 0.8 1.0 

x 

Figure 2: The limit distribution of r]i/K as K — > oo, given by the PDF of 
Beta(y, y) (cf- EID for several values of m. 

Proposition 4.2. In the limit K — > oo we have for all m > 0, 

( -7-, -7-) — >■ (5, 1 — 5) m distribution , (50) 

where B G [0, 1] is a continuous random variable with Beta y) distribu- 
tion and PDF 

f B {x) = ^^x^-\l-xr^\ X6[0,l]. (51) 

Proof. Using ([191) we get as fT — )■ 00 and n/lf — » x G [0,1] for the 
asymptotic form of the numerator of (149 p 

K m - 2 x m/2 - l {l-x) m/2 - 1 . 

For the denominator we get the integral 

K m-2 f 1 m/2-l (l _ )m/2-l g rf = ^-i jH 

io T(m/2) 2 ' 

using the representation £>(r, s) = p^p^) for the Beta function. Thus we 

have that K^(r]i — n) — > /^(a;) converges to the PDF of the Beta(y, y) 
distribution. □ 
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We see that for m < 2 one site contains most of the particles while for 
m > 2 both sites are likely to have around K/2 particles. The boundary 
case is m = 2, where the particles are distributed uniformly among the two 
sites. This is a standard property of the symmetric Beta distribution and 
is illustrated in Figure [2j In the limit m — > we recover the degenerate 
distribution (T48"j) . 

Remark 4.1. a) The result ( f^g| ) can be immediately generalized to a finite 
set An = {1, . . . , N} of N > 2 sites. In the limit m — > we have for 
all K GN 

1 - 

1=1 

where = (.., 0, 1, 0, ..) G ~R N is the standard unit vector in direction i. 

b) In the absence of diffusion for m = the inclusion process has in 
general many absorbing states which exhibit several isolated piles of 
particles. However, if the p(i, j) > for all i, j £ An, then all absorbing 
states have exactly one pile containing all the particles. The stationary 
measures are then all possible mixtures 

N 

2_. a i$Kei with cti E [0,1] and ^^cti = 1 • (53) 

i=l i 

The limit result / T5^j) leads only to the symmetric mixture, due to ho- 
mogeneity and ergodicity of the process for m > 0. 

Connection to zero-range processes. 

This result is slightly different from most previous work on homogeneous 
zero-range condensation, which is mostly discussed in the limit of infinitely 
many particles [13J or the thermodynamic limit [HIE]. In this case, above a 
certain density or particle number all sites have heavy-tailed distributions and 
condensation is a consequence of large deviation properties of such random 
variables, as discussed in detail in [14]. 

For the inclusion process we discuss the two extreme cases of a finite 
and an infinite lattice (see next section), in the limit of a vanishing system 
parameter m — > 0. The distributions of the occupation numbers always have 
exponential tails due to the cut-off (jSj), which disappears in the limit in a 
non-uniform way. This is very similar to results in [TB], where a parameter 
was varied together with the system size in a joint limit. Analogous results 
are phrased here in terms of the law of large numbers in the next section. 
Size-dependent system parameters have also been studied in [19], which can 
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lead to a cut-off similar to (|44[) and a typical maximal cluster size also in 
zero-range processes. 

As a further difference to zero-range condensation, there is no non-trivial 
critical density p c for the distribution of sites outside the maximum in the 
inclusion process. In fact, in the limit m — > all N particles condense on 
a single site, which corresponds to p c = and is an absorbing state for 
the dynamics with m = 0. This is related to results on zero-range processes 
where the jump rates vanish in the limit of infinite occupation number, which 
has been studied in [20J and more recently also in [2T] . 

4.3 The infinite- volume limit 

For finite systems with a fixed number of particles the exponential part of the 
product measures that leads to a cut-off for large n (cf. (114"]) ) did not play any 
role due to cancellation, but will be of importance for infinite systems. For 
simplicity we consider stationary configurations of the symmetric inclusion 
process (SIP) on the infinite lattice A = N which leads to a family of iid 
random variables rji,r]2,--- with distribution i/p (1401) . In this context the 
condensation phenomenon for m — > can be formulated as a breakdown of 
the usual law of large numbers. 

For every m > by definition E(r/i) = p and a usual law of large numbers 
holds, i.e. 

1 K 

Sk := — 7 rji p a.s. as K — > oo . (54) 

8=1 

On the other hand, rji — >■ as m — > in distribution, and therefore we have 
for all K G N even for the unnormalized sums 

K 

rji — > in distr. as m — > . 

i=i 

This implies that the limiting behaviour of the empirical mean as K — > oo 
and m — > depends on the order of limits. Thus we are interested in the 
joint limit K m — > oo as m — > to identify the scale on which the law of 
large numbers changes behaviour. It turns out that there are two interesting 
scales for K m . 

Proposition 4.3. Let K m = — — ■ Then as m — >■ we have (in distr.) 

"* mlogm > ' 

A K m ■= J2 i 1 ~ 6o ^) — > \ W S > K m /K m — )■ 5 G (0, oo) , (55) 

t=l [ OO , K m > K m 
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where W$ ~ Poi(S/2) is a Poisson random variable with mean 5/2. In the 

last case, A Km = + o(l)) . 

Furthermore, on the larger scale l/m> K m we have 



S K m = -rr- ^2 * — ► { X 7 > -> 7 e (0, oo) , (56) 



L m 



K m ( , i<" m < 1/m 

^ % — > l A 7 , K m m -> 7 ( 
»=i ( p , i^m > l/m 

where X y ~ GammaQ, ^) zs a Gamma random variable with mean p. 
Proof. Denote the probability of rji > by 

p w := 1 - „<»0(0) = logm(l + o(l)) = + o(l)) , (57) 

with asymptotics for m — >• 0. Then 1 — 5o jT]i ~ Be(p m ) are i.i.d. Bernoulli 
random variables and therefore A^ m ~ Bi(K m ,p m ) is a Binomial with 

P(A^ m = n)= ^^(l-p m )^-", 

counting the non-zero contributions to the sum Sjc m - p m — > as m — > 
with asymptotics given in (157)) . and (153]) is a well-known scaling result for 
Binomial r.v.s. Since the rescaled random variables (1 — So jT]i )/p m have mean 
1, we have by the ususal law of large numbers 

± ^ml J m - Iv m I'm 

This holds whenever K m p m — > oo or, equivalently, K m 3> Km since the sum 
will have infinitely many non-zero contributions, and implies that Ak = 

Analogous to fl38|) we get for the characteristic function of 77^ 

m/2 



X,(t)=E(e^) = ( T 



For the rescaled sum Sx m of iv" m independent r.v.s we get 
X s(t) = Xv(t/K m f m = [1 + ^ (1 - e^)j 

where we used p = p m (4>) = y yt^ as i n Q55D to fix the density. As K — > 00 
we have for all complex z 7^ 

1 - z 1 '* = -1 log z (l + o(l)) . (58) 
20 



This leads to the asymptotics 



2p 



) 



K m m/2 



Xs(t) 



1 - 



K m m 



it 



(l + o(l)) , 



since the correction terms from (1581) are of order 1/K m <C l/(mK m ). There- 
fore, asm-fO 



which implies the weak law of large numbers in the two extreme cases of (|56|) . 
In the intermediate case mK m — > 7 we have 



This result leads to the following interpretation for the limiting behaviour 
of Sx m as m — t- 0. 

a) K m <C K m : There are no non-zero contributions to Sjc m and even the 
unnormalized sum K m SK m — > 0. 

b) K m ~ K m : There is a finite (Poisson distributed) number of non-zero 
contributions to Sx m , but still Sx m — > 0. Since the law of these con- 
tributions becomes degenerate as m — » (cf. (145]) ) we have no scaling 
law for K m S Km . 

c) K m K m < 1/m: SK m has an infinite number of non-zero contribu- 
tions, but still vanishes as m — >■ 0. 

d) ~ 1/m: S*x m has a random limiting value (Gamma distributed) 
with mean p, and infinitely many non-zero contributions. This inter- 
polates between the deterministic limits and p, as shown in Fig. [3j 

e) K m ^ 1/m: The usual weak law of large numbers holds, i.e. Sic m — > P 
as m — > 0. 

If we interpret rji,rj2, ■ ■ ■ as a configuration of the inclusion process, this 
result gives detailed information about the structure of such configurations 
as m — > 0. They are in direct analogy to results in [TS] on a particular 
zero-range process, which have just been formulated in an inverted fashion 
corresponding to a parameter itlk — > in the limit K — > 00. 




, mK m 
, mK rn — > 00 




which is the characteristic function of a Gamma 



(tt, — ) random variable. □ 
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Figure 3: The limit distribution of Sx m as m — > with mK m — > 7, given by 
the PDF of a Gamma Q, ^) random variable (|56|) . In all cases p = 1, and 
increasing 7 (3 values shown) interpolates between the deterministic limits 
and p. 

5 The Brownian energy process 

In pQ we introduced the Brownian energy process with parameter m > 
(abbreviation BEP(m)), and explained how, for integer values of m it is 
related to the Brownian momentum process with m momenta per site. 

More precisely, the BEP(m) is an interacting diffusion process on Qn = 
[0, oo) 1 '"''^ with generator 

Lf(x) = f> X * +1 (^-^) 

- 2m(x * - {-k - ^b) ' (59) 

where for a configuration of "energies" 1 G Hjv, Xj denotes the energy at site 
(1 V), 

In [3] we introduced an asymmetric version of the Brownian momentum 
process. This model was later studied in [3]. Motivated by this asymmet- 
ric modification of the Brownian momentum process, we now introduce an 
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asymmetric version of BEP(m) via its generator 



L f( x ) = 5^4x^+1 




(60) 



We focus on a one-dimensional nearest-neighbour lattice as for the ASIP 
in Section 3, but the definition could of course be generalized to arbitrary 
geometries. Obviously, the total energy f(x) = YliLi x i ls conserved, and for 
E > the process has a drift to the left, which can most easily be seen from 
the stationary measures discussed in the next section. 

5.1 Condensation in the ABEP 

We first consider m = 2, E > 0, and two sites. This is the simplest case 
because the marginals of the stationary distribution are exponential, which 
makes explicit computations simple. The generalization to m > and more 
sites is easy. 

The generator, written in the variables (xi,^) = : (u,v), then reads: 



The adjoint (in L 2 (R, dx)) is given by (the closure of the operator) 



L 




L* 




As an ansatz for the density of the stationary distribution we put 

f(u,v) = abe~ au e- bv 



(61) 
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with a, b > 0. Plugging this in the equation for the stationary density L*f = 
gives 

Auv(b - a) 2 + 4(v - u)(b - a) + 2E(u - v) - 2£w(6 - a) = , 
which leads to 

6 = a+| (62) 

and 

f(u, v) = a(a + E/2)e- au e- av e- Ev/2 . (63) 

In order to state our condensation result, denote by (Uk, Vk) the pair (U, V) 
with probability density (1631) conditioned on U + V = K. We then have the 
following result, which should be thought of as the analogue of Theorem 13. 11 
but now in continuous state space setting. 

Theorem 5.1. a) As K — )■ oo, Vk converges in distribution to a random 
variable with exponential distribution with parameter E/2, i.e., with 
probability density (E/2)e~ u( - E ^ . 

b) As K —7- oo ; Uk/K — > 1 almost surely. 

Proof. The proof is a direct computation. Put A = a + E/2, X' = a, then 
A > A', A — A' = E/2. First note that the distribution of U + V has probabily 
density ^(e"^ - e~ Xx ). 

Next, the conditional density of V given U + V — K is given by 

AA'(A - \>) e -^ e -\>{K-u) ^ (A - AQe-( A - A > 
\X'(e- x ' K - e- XK ) ~ i _ e -(x-\')K 

which converges, as K — > oo to (A — \')e~( x ~ x '^ u , implying statement a) of 
the theorem. To prove statement b): choose < 8 < 1, then 



F(U < (1-6)K\U + V = K) 



&- 6)K XX' {X - X')e-* x 'e-( K -^ x dx 



XX'(e~ x ' K - e~ XK ) 
(A - A') J^ 5)K e-< x ~ x 'Mx e (A-A')^d-5) _ i 



e {\-\')K _ l e (\~\')K _ I 

K l_ e (A'-A)K 
\)K 



as K — > oo. As in the proof of Theorem 13.11 the bound is summable in K 
if we choose 5 = 1/yfK and Uk/K < 1 by definition, which implies almost 
sure convergence. □ 
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To generalize the previous computation to the case of N sites and general 
parameter m > 0, it is easy to check along the lines of the proof of Theorem 



with 



12.11 that the process with generator (|60|) has a stationary measure which is 
a product of Gamma distributions with identical shape parameter m and 
site-dependent location parameter. More precisely, the PDF is given by 

N mil m/2-1 _„. T . 

/(*„-,*„) = n ' r' (m /2) (64) 

(i-l)E . . 

a i = a+ 1 ; (65) 

for i G {1, . . . , N}. After conditioning on the sum X\ + . . . + Xx = K 
we find, again by simple explicit computation, in the limit K — > oo that 
X\/K converges to 1 almost surely, and that for i = 2, . . . , N, the law of Xi 
converges to a shifted Gamma distribution with density 

lim f Xt \x 1+ ...x N =K(xi) = Ctx? e ? (66) 

where Cj = (^e)^"/r(y) is a normalization constant. 

The interpretation of this result is the same as in the discrete case for 
the ASIP. Here almost all energy concentrates on the lattice site with the 
heaviest tail in the stationary distribution. 

5.2 Generalizations 

Exactly as in the case of condensation in the ASIP (section 3.2), we can for- 
mulate a more general condensation result for independent random variables 
Xi, . . . , Xn with values in [0, oo) and marginal densities 

fx 1 (x) = ^- ) e- x ^w l (x)e^ . (67) 

where < Ai < min^ 2 \j. Here a notation with so-called chemical potentials 
/iGlis more convenient than the fugacity variable = e M used for the SIP, 
and values — oo < \l < fi c := Ai are possible. The normalization 

POO 

Zi(n) — / e~ XiX Wi(x) e^ x dx 



o 



is finite for /i < fi c , and for indices i < N also Zi(fj, c ) < oo. The Wi : [0, oo) — > 
[0, oo) are subexponential in the sense that for all y G K 

lim W{X + V) = 1 . (68) 



25 



The proof of this result follows the same steps as the proof of the analogous 
discrete result, except that we have to replace sums by integrals. As this is 
a straightforward extension, we leave the proof to the reader. 



Theorem 5.2. Denote by (Y 1 , . . . , Y^) the random variables . . . , X 
conditioned on Xi + . . . + -Xjy = K . Then under the above conditions §6% 
and we have as K — >■ oo: 

a) Condensation on the site with the heaviest tail, i.e. 

Y K 

— > 1 almost surely ; 

b) Convergence to the critical distribution with fi = fi c for other sites, i.e. 

(F 2 X , . . . , Yx) — > (Y 2 , . . . , Yjv) in distribution , 
where the Y^ are independent with densities 

fvM = -^—e-^e^ Wl (y) . 

Remark 5.1. In the limit m — > 0, also for spatially homogeneous Brownian 
energy processes there will be a condensation phenomenon as m — > com- 
pletely analogous to the results in Section 4 for the inclusion process. Indeed, 
for a fixed average energy p > (taking ai = m/(2p) in the marginal 

densities of the stationary product measure are 



(\ m/2 
— ] x m / 2 - 1 e -mx i /(2p) 



Analogous to Theorem J^l one can easily show that this implies 

nx t <5)= [ fxXx^dx^l 
Jo 

for all 5 > as m — )■ , so that Xi — > in probability. Further, all statements 
following from Theorem 4-1 in Section 4 can be derived in an appropriate 
version for continuous variables. 
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6 Conclusion 



We have studied condensation phenomena for random variables with expo- 
nential tails, which arise in the inclusion process and related particle systems. 
In general, condensation can be due to the presence of subexponential tails 
resulting from a strong particle attraction, which has been studied in detail 
in the context of zero-range processes [SJ [TT], [T21 [321 [HJ E3- For exponential 
tails considered in this work, the attraction between particles alone is not 
strong enough and a second ingredient is needed for condensation. 

One possibility are spatial inhomogeneities, which will lead to a non-zero 
fraction of the particles to cluster on the sites with the heaviest tails in the 
limit of infinitely many particles. Our result on this in Section 3 applies in 
great generality, extending also previous related work on zero-range process 
[3 El El HO]. For homogeneous systems, varying a system parameter can in- 
duce condensation for fixed total particle density as studied in Section 4 for 
the inclusion process. Previous results in that direction include [TBI EE] for 
zero-range processes and also [22J for a continuous mass model. The Brow- 
nian energy process studied in Section 5 provides an interesting example 
where both versions of condensation can be studied in a system with contin- 
uous state space and dynamics. Condensation for continuous variables has 
been studied before in the random average process [22] and mass transport 
models [221 121]) all of which use a discontinuous redistribution of mass (or 
energy) following a jump process. 

To summarize, inclusion processes and related systems such as the BEP 
provide a rich class of models that exhibit condensation phenomena of several 
kinds in the presence of exponential tails, the description of which applies also 
in more general situations. For inhomogeneous models we have focused on 
finite systems, and a further question would be to consider thermodynamic 
limits where, for example, inhomogeneity is due to random disorder as studied 
in [Zl El El HO] for zero-range processes. In the homogeneous case it would be 
of great interest to exploit duality in the SIP and BEP to get results on the 
dynamics of condensation. 
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